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DYNAMICS  OF  SPACE  CHARGE  WAVES  IN  THE 
LASER  BEAT  WAVE  ACCELERATOR 

I.  Introduction 

The  laser  beat  wave  accelerator  is  one  of  a  number  of  laser  driven 

1  2 

particle  accelerating  schemes,  *  which  is  currently  receiving  considerable 
attention.  Existing  types  of  synchrotron  and  linear  accelerators  are  nearing 
their  economic  limits  in  going  much  beyond  a  few  TeV  in  energies.  The 
availability  of  high  power  laser  beams  (£  10  watts)  with  electric  fields  as 
high  as  10^V/cm  brings  about  the  possibility  of  using  these  high  fields  to 
accelerate  particles.  Direct  use  of  these  fields  for  continuous  particle 
acceleration  is  not  possible  due  to  the  transverse  polarization  and  rapid 
oscillation  of  the  fields.  A  number  of  laser  driven  acceleration  schemes  have 
been  suggested  over  the  past  two  dozen  years. 

.  Before  describing  the  laser  beat  wave  accelerator  concept  a  brief 
description  of  some  of  the  other  generic  laser  acceleration  concepts  are 
mentioned.  These  include  the  inverse  free  electron  laser  accelerator,  the 
grating  accelerator,  the  inverse  Cerenkov  accelerator,  the  cyclotron  resonance 

accelerator  and  the  use  of  high  gradient  short  wavelength  structures. 

3—5 

In  the  inverse  free  electron  laser  accelerator  scheme,  an  electron 
beam  together  with  an  intense  laser  pulse  is  propagated  through  a  spatially 
periodic  magnetic  field  known  as  a  wiggler  field.  The  wiggler  period  and 
laser  wavelength  are  such  that  the  transverse  particle  velocity  due  to  the 
wiggler  field  is  in  phase  with  the  transverse  electric  field  of  the  laser 
radiation.  By  appropriately  contouring  both  the  wiggler  amplitude  and  period, 
the  injected  particles  can  be  continually  accelerated.  The  inverse  of  this 
process  has  been  used  to  generate  radiation  and  is  the  well-known  free 
electron  laser  mechanism. 

The  grating  accelerator  mechanism6  relies  on  the  fact  that  when  electro¬ 
magnetic  radiation  propagates  along  a  diffraction  grating  a  slow  electro- 
Manuacript  appro  red  September  10,  1984. 
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magnetic  surface  wave  is  excited  along  the  grating's  surface.  This  scheme 
utilizes  the  slow,  phase  velocity  less  than  the  speed  of  light,  electro¬ 
magnetic  wave  to  trap  and  accelerate  a  beam  of  injected  electrons. 

7  ft 

The  Inverse  Cerenkov  accelerator  ’  approach  takes  advantage  of  the  fact 
that  the  index  of  refraction  of  a  neutral  gas  is  slightly  greater  than 
unity.  The  laser  radiation  within  the  gas  has  a  phase  velocity  less  than  the 
speed  of  light  making  it  possible  to  trap  and  accelerate  an  injected  beam  of 
particles. 

9  10 

In  the  cyclotron  resonance  accelerator  mechanism  *  an  electron  beam  is 
injected  along  a  uniform  magnetic  field  together  with  a  parallel  propagating 
iaser  beam.  Because  of  a  self  resonance  effect,  the  phase  of  the  electron's 
transverse  velocity  can  be  synchronized  with  the  radiation  electric  field. 

This  synchronism  is  maintained  throughout  the  acceleration  length. 

The  high  gradient  short  wavelength  structure  concept11  is  basically  a 
scaled  down  version  of  a  conventional  slow  wave  accelerator  structure. 
Radiation  power  sources  in  the  centimeter  wavelength  range  appear  appropriate 
for  this  approach.  The  potential  advantage  of  this  scheme  is  that  due  to  the 
short  wavelength  employed,  relatively  low  radiation  energy  per  unit  length  is 
needed  to  fill  the  small  structure,  and  breakdown  field  limits  appear  to  be 
higher. 

The  laser  beat  wave  accelerator  concept  is  a  collective  acceleration 
scheme  which  utilizes  a  large  amplitude  plasma  wave  with  phase  velocity 
slightly  less  than  the  velocity  of  light  to  accelerate  charged  particles.  The 
large  amplitude  plasma  wave  is  generated  by  the  nonlinear  coupling  of  two 
intense  laser  beams  propagating  through  the  plasma.  *  In  this  process 

the  two  laser  beams  with  frequencies  oi^,  and  corresponding  wavenumbers  kj, 

k2  couple  through  the  plasma  to  produce  a  ponderoraotive  wave  with  frequency 


ai,  -  oi-  and  wavenumber  k,  -  k«.  See  Fig.  1.  If  oi,  -  oi_  *  oi  ,  the  plasma  wave 
I  i  *  £  I  2  p 

will  initially  grow  linearly  in  time.  If  the  laser  frequencies  are  much 
greater  than  the  ambient  plasma  frequency  oi^  then  the  phase  velocity  of  the 
ponderomotive  wave  is  nearly  equal  to  the  group  velocity  of  the  laser  wave. 

In  this  scheme  a  beam  of  injected  electrons  with  axial  velocity  close  to  the 
plasma  wave  phase  velocity  can  be  accelerated  until  synchronism  is  lost. 

A  potentially  attractive  variation  of  the  plasma  beat  wave  accelerator  is 
the  surfatron  scheme.  In  the  surfatron  configuration  a  transverse  magnetic 
field  is  externally  applied  permitting  the  accelerated  particles  to 
effectively  E  x  B  drift  in  a  direction  transverse  to  the  laser  propagation 
direction.  In  this  configuration  the  electrons  can  remain  in  phase  with  the 
plasma  wave  allowing,  in  principle,  higher  electron  energies  to  be  achieved. 

In  this  paper,  we  analyze  the  build-up  of  the  plasma  waves,  extending  the 
analysis  of  Rosenbluth  and  Liu.^’^  Our  nonlinear,  fully  relativistic 
treatment  of  the  plasma  wave  includes  mismatching  of  the  laser  beat 
frequencies  to  the  plasma  frequency,  applied  transverse  magnetic  field  as  well 
as  time-dependent  laser  pulses.  The  resultant  equations  in  Section  II 
describe  the  growth  of  the  excited  plasma  waves  up  until  saturation.  We  find 
that  the  effective  plasma  frequency  is  increased  from  the  ambient  plasma 
frequency  when  a  transverse  magnetic  field  is  applied.  On  the  other  hand,  the 
effective  plasma  frequency  can  be  decreased  as  the  transverse  motion,  induced 
by  intense  laser  powers,  becomes  relativistic. 

Making  the  weak  laser  power  approximation  in  Section  III,  we  obtain 
analytical  results  for  the  initial  growth  rate  for  the  large-amplitude  plasma 
wave,  the  maximum  accelerating  electric  field,  the  laser  beat  frequency 
requirement  and  the  corresponding  phase  velocity.  As  the  combined  laser 
powers  (measured  by  e)  increases,  the  maximum  amplitude  of  the  olasraa  electric 


field  increases  as  ae  .  The  variable  a  is  a  function  of  frequency  mismatch 
between  the  laser  beat  frequency  and  the  effective  plasma  frequency.  The 
relativistic  effect  associated  with  the  transverse  motion  is  incorporated  in 
the  parameter  e.  In  the  limit  of  exact  resonance  and  nonrelativistic  motion 

*s  / 

in  the  transverse  direction,  we  confirm  the  results  of  Rosenbluth  and  Liu. 

For  the  purpose  of  accelerating  electrons,  it  is  desirable  not  only  to  have  the 

largest  accelerating  electric  field  but  also  the  phase  velocity,  v  ,  ,  less 

ph 

than  but  close  to  the  speed  of  light.  We  find,  however,  that  as  the  amplitude 

of  the  plasma  wave  gets  larger,  the  corresponding  phase  velocities  become 

2/3 

slower  i.e.,  v  -  v  scales  as  e  .  where  v„  is  the  group  velocity  of  the 
g  ph  g 

lasers  in  the  plasma.  The  maximum  accelerating  field  for  a  given  laser  power 

parameter  e  is  achieved  when  the  laser  beat  frequency  is  less  than  the 

effective  plasma  frequency.  The  difference  between  the  effective  plasma 

2/3 

frequency  and  the  optimal  laser  beat  frequency  is  proportional  to  e  .  The 

-2/3 

plasma  build-up  time  is  proportional  to  e  .  The  analytical  results  also 
shows  that  the  transverse  magnetic  field  has  little  effect  on  the  plasma 
dynamics  to  the  lowest  order. 

In  Section  IV,  we  obtained  numerical  results  from  the  full  treatment  and 
compared  them  with  the  analytical  results.  In  the  weak  laser  power  limit,  the 
results  are  in  excellent  agreement.  We  observe  the  wave  steepening  phenomenon 
caused  by  relativistic  effects.  As  the  laser  power  increases,  the  waves 
eventually  break  and  become  turbulent.  We  show  that  for  laser  pulse  lengths 
much  longer  than  the  plasma  wave  build-up  time,  the  amplitude  and  the  phase 
velocity  of  the  plasma  waves  are  oscillatory.  By  applying  lasers  with  pulse 
lengths  approximately  equal  to  the  plasma  wave  build-up  time,  both  the 
amplitude  and  phase  velocity  can  be  approximately  maintained  at  a  constant 
value. 
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II.  Plasma  Beat  Wave  Accelerator 


Our  model  consists  of  a  spatially  one-dimensional  plasma  containing 
infinitely  massive  ions.  Initially  the  plasma  is  assumed  to  be  cold,  uniform 
in  density  and  stationary.  Large  amplitude  plasma  waves  are  excited  by  the 
ponderomotive  force  associated  with  the  two  laser  pulses.  Using  a  Lagrangian 
formulation,  the  temporal  evolution  of  the  plasma  wave  over  a  single  spatial 
period  is  studied  at  a  fixed  axial  position. 

This  analysis  treats  the  following  topics:  1)  nonlinear  behavior  of 
plasma  waves,  2)  relativistic  effects,  3)  effect  of  finite  duration  laser 
pulses,  4)  mismatching  of  the  laser  beat  frequency  to  the  plasma  frequency, 
and  5)  the  effect  of  an  applied  transverse  magnetic  field. 

A.  Non-Linear,  Relativistic  Plasma  Dynamics 

The  vector  potential  associated  with  the  linearly  polarized  laser 
pulses  within  the  plasma  is 

A 

A  (  z ,  t )  *  Z  A,(z  -  v  t)  cos(k.z  -  oj.t  +  <i> .  )e  ,  (1) 

1-1,2  1  g  i  i  i  x 

where  A^(z  -  v^  t)  denotes  the  pulse  amplitude  of  the  itll  laser,  »  oip,  o>^ 
is  the  laser  frequency,  mp  -  (4ir|e|  n^/n^)  is  the  ambient  plasma  frequency, 
nQ  is  the  ambient  electron  density,  and  v^  -  (u>^  -  -  k^)  is  the  group 

velocity  of  the  laser  radiation.  In  our  model,  we  assume  k^  to  be  constant 
and  A^  specified,  i.e.,  the  imposed  laser  fields  are  assumed  to  be  unaffected 
by  the  plasma  density  modulations.  This  is  a  good  approximation  if  »  o)p, 
je^/fn^c  )  <  1  and  5n/nQ  <  1*  where  5n  is  the  plasma  density  modulation. 

Also  included  in  our  model  is  an  applied  transverse  static  magnetic  field 
which  is  represented  by  the  vector  potential 
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A  *  B  z  e  .  (2) 

*"*o  o  y 


The  electric  and  magnetic  fields  in  terms  of  the  vector  potential  are 
1  3 

£  *  ~  —  Al  and  3  -  7  *  (A^  +  Aq)  respectively.  In  our  configuration,  the 
electric  field  vector  is  parallel  to  the  transverse  magnetic  field,  E^SB^. 

In  our  one-dimensional  model,  the  transverse  particle  dynamics 
possess  two  constants  of  motion,  i.e.,  the  canonical  momentum  in  the  x  and  y 
direction. 


*1  •  -  ¥<41.  ♦  V  *  Si 


(3) 


Assuming  that  *  0  prior  to  the  arrival  of  the  laser  pulses,  it  follows  that 
the  electron's  transverse  momenta  are  given  by 


p  -  M  at  •  e 
Kx  c  > 


(4) 


and 


B(z-z) , 
CO  o 


where  zQ  is  the  initial  axial  position  of  the  electron. 


(5) 


The  momentum  equation  in  the  axial  direction  can  now  be  written  as 


dt 


l  |  3A  3A 

-lei  E - Lai  (p  _JL  +  P  3-2.), 


z  Y®0c  x  9z 


'  y  3z 


(6) 


where  &z  is  the  self-consistent  plasma  induced  electric  field,  Ax  “  *  ex 

and  Ay  ■  BQz.  The  particle  energy  equation  is 


*4* 


(7) 


,  .  ,  3A  3A 

il  =  -LSI-  (P  _2L+  P  — Z.)  _  _  e| 
dt  23  ^  x  3t  y  3t  ;  22  z  *# 

1®  c  yh  c 

o  o 

Combining  (6)  and  (7),  we  obtain  the  equation  governing  the  axial  motion  of 
the  electrons. 


(l  .  3v/2 
|ej_  (1  V 

V  (1  +  U2  +  u2)1'2 
x  y 


c  o  -  B  )  8  B  8  2 

- 5 —  (z —  +  _£  L_)  ry2 

2(1  +  U2  +  U2)  3z  c  3t  x 

x  y 


where  0z  =  v^/c,  Ux  -  ^  ^  *  Vand  U  -  -  B  (*-*  ). 

me  7  m  c 

o  o 


Neglecting  nonresonant  terms  we  find  that 


Ux  +  Uy  =  ^z~2o^2  +T  ^al  +  a2?  +  ala2  C0s  ^Akz  ~  Aut  +  (9) 

c 


where  n  =»  |e|B  /m  c  is  the  cyclotron  frequency,  a=|e|A, (z-v  t)/m  c  , 
o  o  i  1  g  o 

Ak  *  k^  -  k^,  A<d  *  oi^  -  oj^,  A$  ■  and  8^  *  Aw/cAk.  The  difference  in 

the  laser  frequencies  is  assumed  to  be  close  to  the  plasma  frequency, 

Am  ■»  ujp.  Substituting  (9)  into  (8),  yields 

dB*  _  ]e[  <>-B2)  a2  . 

dt  me  y  z  2  c  o 

o  l 


0-6  )  6  a  17  9 
+  c  ~  2  %  +  r  ar^2(ai  +  aP  +  ai  a2  cos(Akz 

iyL 


-  Acot  +  Aifc  )  j  , 


where 


y.  «  (1  +  U2  +  U2)1/2 
1  x  y 


7 


2 

*  11  +  — j  (z-ZQ)2  +  (a2  +  a22)  +  a^a2  cos(Akz  -  Acot  +  A4>)]  ^2 

c 

is  a  measure  of  the  magnitude  of  the  transverse  oscillations  induced  by  the 
iaser  beams. 

B.  Transformation  to  Lagrangian  Variables 

Equation  (10)  is  expressed  in  terras  of  Eulerian  independent  variables 

z  and  t.  It  proves  convenient  at  this  point  to  perform  a  transformation  to 

23  24 

Lagrangian  variables,  *  because  the  plasma  induced  field  Ez  takes  on  a 
particularly  simple  form.  In  the  Lagrangian  variables,  the  independent 
variables  are  zQ  and  t,  such  that 

t  -  t,  (12) 


z  -  zq  +  5(zq,t), 

where  zQ  is  the  initial  position  of  the  electrons  at  x  ■  0  and  E(zq,t)  is  the 
axial  displacement  at  time  T  relative  to  its  initial  position  zQ.  The  axial 
electric  field  written  in  terras  of  the  new  variables  (zq,t)  is 

Ez(zo,T)  -  Air  |  e  |  nQ  C(zq,t).  (13) 

The  equation  governing  the  plasma  oscillations  can  now  be  written  completely 
in  terms  of  S(zq,t),  i.e.. 


8 


3a 


V. 


U 


a 

a 


3a 

+  (a^  +  a^  cos  Ai|»)  — 

I 

-  (3  l-  |/c)  cAuj  a1  a2  sin  A<j>f  , 


where  I  *  3£(zq,t)/3t,  a  *  zq  +  £  -  v^t  and  AiKzq,t)  =  Akf  zQ  +  £(zq,t)I 
-  Acot  +  A<fr. 

The  imposed  laser  field  and  the  applied  transverse  magnetic  field 
modify  the  plasma  frequency  and  hence  we  define  an  effective  plasma  frequency 


where 


'io 


12  2 
[1  +  ,(a.  +  a_  )] 

2  1 ,max  2, max 
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(16) 


The  effective  plasma  frequency,  '0g,  is  the  relativistic  upper  hybrid 
frequency.  An  increase  in  the  transverse  oscillations  results  in  a  decrease 
in  the  velocity  of  the  axial  oscillation,  which  in  turn  leads  to  an  effective 
reduction  in  the  plasma  frequency.  On  the  other  hand,  the  transverse  magnetic 
field  results  in  an  increase  in  the  effective  plasma  frequency. 

C.  Normalized  Plasma  Wave  Equation 

To  further  study  the  dynamics  of  the  plasma  oscillations,  we 
normalize  the  parameters  in  the  following  manner, 


T  -  Aojt, 


(17) 


.  M  E 

2  m  z 
U)  o 

P 


Ak!i, 


9 


Zn  =  Akz^ , 
o  o 


f  =  Aoj/o)  , 
s 


G  =  (^/Yi02)  («p2/Tio)_1. 


The  evolution  of  the  plasma  oscillation  is  now  given  by 


+  - - I -  (i  -  g2  g2)  Ii°  [(l  -  g2  g2)1/2  +  g  E 

_Zv.  _  *  g  Y  c  vJ 


f  (1  +  G) 


(1  -  82  g2)  (l  -  z1  g)  9a. (a) 

- S-j-j - lUl  +  a2  C0S 


2i  t 


3a2(a) 

+  (a2  +  a i  cos  A\|>)  — ^ - al  a2  sin  A,H  » 


where 


11  +  72“  TTg  Y1o  e2  +  i  (al  +  d2}  +  ala2  COS  A*i  1/2’ 


Ai|)  *  E  -  T  +  ZQ  +  A<J>, 


a  *  Zq  +  E  -  T  and  E  »  9E(Zo,T)/9T.  The  total  relativistic  mass  factor 


associated  with  the  plasma  motion  becomes 


-2  i 2  ~l/2 

T(Zo>T)  -  Y±  (1  -  8  Z  £2) 


’*  *»  ••  '•  *.  N  .  .  *•  ; 
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Since  the  plasma  oscillation  amplitude  Is  single  valued,  our  model  is 
valid  prior  to  electron  trapping,  i.e.,  trajectory  crossing.  Trajectory 
crossing  occurs  when  the  amplitude  | E I  is  approximately  unity. 

The  equation  of  motion  described  by  Eq .  (18)  is  fully  relativistic  in 
both  the  axial  and  transverse  directions  and  permits  us  to  model  laser  pulses 
with  a  beat  frequency  other  than  the  plasma  frequency.  Also  by  varying  the 
parameter  G  associated  with  the  transverse  magnetic  field  we  can  also  model 
the  surfatron  configuration. 


11 


III.  Analytical  Theory  in  the  Weak  Laser  Field 

In  this  section  we  derive  the  scaling  laws  appropriate  for  weak  laser 
fields,  i.e.,  electron  transverse  motion  is  mildly  relativistic.  Here,  we 
consider  the  case  where  the  two  laser  powers  are  constant  in  time  and  define 
the  small  parameter 


£ 


a  a 
12 


1  +  I*al2  +  a2  ^ 


(22) 


where  the  denominator  of  e,  is  associated  with  the  relativistic  motion  in  the 

2 

transverse  direction,  Y1q.  When  aj  *  a2»  e  Is  proportional  to  the  laser 
powers. 

Assuming  E  to  be  a  slowly  varying  function  of  time,  it  can  be  represented 
in  the  form 


E(Zq,T)  -  AE(T)  sin  (Zq  -  T  +  0(T)  +  A<fr) , 


(23) 


where  AE(T)  and  9(T)  are  slowly  varying  functions  of  time  compared  to  the 
plasma  oscillation  period.  Furthermore,  we  assume  that  u>  ,  w  »  a>  ,  so 

12  p 

that  Aio/cAk  =  1. 

A.  Small  Parameter  Expansion 


I 


Expanding  Eq.  (18)  in  terms  of  the  small  parameters  e  and  E 
obtain  the  following  equation  for  E, 


we 


E 


|  sin  A*  +  F 


(24) 


where 


I 

f(zq  »t) 


“  ^2^^  0  “  y  S^)  (!“■§■  cosAil/)  -  (1  “  e  cos  Aty)]E 


12 


-  y  [E  +  £  cos  A^jsin  &<t>, 


A<l»  *  AE(T)  sin  $(Z  ,T)  -  T  +  Z  +  A4>,  and  t  =»  Z  -  T  +  Ad>  +  6(T).  Using  the 


identities , 


sin  Aiji  -  l  J  (AE)  sinfU  +  1)*  -  6), 

Km— «°  X 


cos  A4<  »  imlm  J^(AE)  cosf  (£+!)#  “  0|, 


(25a) 


(25b) 


and  keeping  only  terms  that  are  proportional  to  sin  $  and  cos  # ,  we  obtain 


( 1  ■  ^7  "  2  af }  AE  sin  $  +  2  ^  cos 


where 


j  Jq(AE)  sin($  -  fi)  -  F,  (26) 


3  2 

F  -  {|  |  AE  +  ij  j^ae)  cos  0  +  J^aE)  cos  20 f  sin 

2 

+  {-^  [(1  2  3G)  +  lj  CAE)  sin  0  -  j-  Jj (2AE )  sin  20|cos  *. 


Using  the  small  parameter  expansions  for  the  Bessel  functions,  i.e., 
JQ(x)  *  1  -  x^/4  and  J^(x)  •  x/2,  Eq .  (26)  becomes 


2  ^  AE  sin  t  -  2  cos  *  ■  AE(i  -  -^)  sin  ® 


c  AE 

+  j  (1  - %  )[sin  $  cos  0  -  cos  $  sin  0j  +  F, 


where 


3  2  2 

F  "  [t  ^2  "  I  (2  ‘  9G)  C0S  9  +  T~  *E  cos(20)]sin  * 


+  It  (2  -  3G)  sin  9  -  7-  AE  sin(20)]cos  ft. 

o  4 


Separating  the  terms  proportional  to  sin  $  and  cos  ft,  we  obtain  two 
simultaneous  coupled  equations  for  the  amplitude  and  phase  variation  of  the 
excited  plasma  wave, 


dAE 

TT 


•li-T&i  +7-^)Jsln  0  +£^sin  20> 


(28a) 


AE 


d9  AE  f2-  1  .  3  AE3  .  .,1  AE2,2  .  ,  9G, ,  „  .  e2 

dT  *  ~2  72  +  I?  ~2  +  ^4  ~  I^“2  +  1  -~2^\^s  0  +  g-AE  cos  20. 

r  r  f  r 


(28b) 


B.  Constant  of  Motion 

Analytical  expressions  for  AE  and  phase  velocity  Vp^  of  the  plasma 

wave  can  be  obtained  by  assuming  1  >>  AE  »  e .  It  will  be  shown  later  that 

2 

this  assumption  is  well  founded.  Neglecting  terms  proportional  to  cAE  and 
2 

c  AE  in  Eq.  (28)  we  obtain. 


dAE 

TT 


■y  sin  0» 


(29a) 


f2  - 

~7 


♦* 


cos  0  + 


3  AE' 

TG 


(29b) 


Multiplying  Eq.  (29a)  by  cos  9  and  Eq .  (29b)  by  sin  0  and  adding  the  two 
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equations  yield 


sin  0,  (32a) 

and 

d9  _ 1  ,,2  9  tV2. 

dT  “  .  ,2  (f  1  +  16  &E  (32b) 

4  f 


where  the  initial  conditions  are,  AE  •  0,  and  8  ■  ir/2.  An  alternate  way  to 

solve  AE  and  0  in  time  is  outlined  in  Appendix  A.  At  exact  resonance  f  *  1  and 

2 

nonrelativistic  motion  in  the  transverse  direction,  i.e. ,  v.  =1  and  c  =  a,a„, 

lo  12 

the  expressions  (32a, b)  agree  with  that  of  Rosenbluth  and  Liu.^ 

C.  Analytical  Results 

We  obtain  analytical  results  for  the  start-up  as  well  as  the  saturated 
regimes  of  the  plasma  wave.  The  plasma  wave  initially  grows  linearly  in  time, 
i.e.,  the  initial  amplitude  of  Ez  in  c.g.s.  units  is  pr^oortional  to  Tt,  where 


r 


6„  ra«c  o 

8  o  2 

7*  -  60 

4  |  e  |  p 


e. 


(33) 
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The  most  interesting  results  are  associated  with  the  saturated  regime.  In  the 
remainder  of  this  section,  we  obtain  maximum  accelerating  electric  fields,  the 
appropriate  laser  beat  frequencies,  the  corresponding  phase  velocities  and  the 
plasma  wave  profiles. 

The  amplitude  of  the  electric  field  is  proportional  to  a  real  root  of 
the  cubic  polynomial  (31).  The  roots  are: 


AE  -  A  +  B, 


A  +  B  ,  A  -  B  r  ~^i 
2  2  ! 


A  +  B  A  -  B  /tt  ,  /  \ 
— - -  / -j  ,  (34) 


where  A  ■  (-  b/2  +  h)1/3,  B  *  (-  b/2  -  h)1/3,  h  -  (b2/4  +  a3/27)l/2, 

2  2 

a  =*  16/3  (f  -  1)  and  b  *  16/3  f  e  cos  9.  For  h  >  0,  there  is  one  real  root 
given  by  AE  »  A  +  B.  When  h  <  0,  there  are  three  real  roots,  and  numerical 
results  show  that  the  relevant  AE  corresponds  to  the  smallest  real  root.  At 
h  =>  0,  AE  undergoes  a  discontinuity. 

For  a  given  laser  power  parameter  e,  the  maximum  AE  occurs  at  h  *  0 
and  cos  0  -  -  1,  giving 


AE 


max 


4  (f) 


1/3 


(35) 


where  e  is  defined  in  Eq .  (22).  The  actual  maximum  field  in  c.g.s.  units  is 


z  ,max 


„2 


m  c 
o 


-e. !/3 

g  "R  (7’  ' 


(36) 


As  the  laser  power  is  increased,  relativistic  effects  on  the  electron 
motion  become  significant.  These  relativistic  effects  cause  the  accelerating 
field  to  maximize  at  a  laser  beat  frequency  which  Is  less  than  the  effective 
plasma  frequency.  For  example  when  the  beat  frequency  is  exactly  equal  to  the 


t 
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effective  plasma  frequency,  i.e.,  f  *  1,  the  maximum  normalized  field  is 


AE  -  (-^-  £)1/3  <  AE  .  (37) 

3  max 

It  can  be  shown  that  the  electric  field  maximizes  to  the  value  in  Eq . 
(35),  when  the  normalized  laser  beat  frequency  is 

V  • 1  - 1  <!  e>2/3-  «»> 

which  corresponds  to  a  laser  beat  frequency  in  c.g.s.  units  given  by 

Aa)opt  *  V1  “  0,54  s2i/3)*  (39) 

During  the  acceleration  process  the  injected  electrons  must  be  nearly 
synchronized  with  the  phase  velocity  of  the  plasma  wave.  To  obtain  the  phase 
velocity,  we  follow  a  null  in  the  plasma  wave.  We  find  that  the  phase 
velocity  of  the  plasma  wave  and  the  associated  relativistic  mass  factor  are 


vg(i  -  d9/dT) , 


(40a) 


and 


V  "  V1  +  2  yl  d0/dT)”1/2’  (40b) 

2  “1/2 

where  d9/dT  is  given  by  (32b)  and  Yg  ■  ( 1  ~  (v^/c)  1  .  As  the  amplitude 

of  the  plasma  wave  becomes  larger,  the  phase  velocity  of  the  plasma  wave 
decreases.  The  phase  velocity  of  the  plasma  wave  at  the  maximum  electric 


field  amplitude  is  a  minimum  and  given  by 


•  »8  <  1  -  l-M  e2/3>- 


(41) 


The  time  it  takes  the  amplitude  of  the  plasma  wave  to  reach  the  first  peak  is 
called  the  plasma  build-up  time  and  in  normalized  units  is  given  by 


T(f,e) 


4  j 


_ de _ 

w/2  f2  -  1  +  (9/16)^’ 


(42) 


where  AE  is  expressed  in  terms  of  0  in  Eq.  (34).  The  equation  above  can  be 

integrated  at  exact  resonance,  i.e.,  f  *  1.  Using  the  fact  that  AE  » 

,  16  _  -.1/3 

(-  y-  e  coso)  at  f  =■  l,  we  can  write 


T(l,e) 


64 

T 


(-3)2/3 

M6J 


e~2/3j 


d0 


k/2  (-  cose) 


173* 


(43) 


Evaluating  the  integral  in  (43),  we  obtain 


?(l,e) 


32  ,-W 3 


8.48c 


-2/3 


(44) 


Taking  e  *  0.01  as  an  example,  the  plasma  build-up  time  at  exact  resonance  is 
29.1  plasma  periods,  i.e.,  T(l,0.01)  ■  29 . 1  (  2tt  ) .  For  e  *  0.1,  the  plasma 
build-up  time  is  reduced  to  only  6.3  plasma  periods.  The  plasma  build-up  time 
shortens  as  laser  power  parameter  c  increases. 

Finally,  we  discuss  the  effect  of  unequal  applied  laser  pulse  powers  on 
the  plasma  wave  dynamics.  Suppose  a  given  value  of  e  is  desired,  the  total 
power  required  for  the  laser  is  at  a  minimum  when  a.  ■  a,  ■  a  where 


The  expression  for  &2  given  and  e  is 


a 


1/2 


(2  +  a^)) 


(45) 
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A  plot  of  a^  +  a^  is  given  in  Fig.  2  for  e  ■  0.1.  The  use  of  unequal  laser 
powers  can  be  advantageous  in  controlling  the  value  of  the  parameter  e. 


IV.  Numerical  Results 


Numerical  examples  are  given  for  the  complete  nonlinear  and  fully 
relativistic  equations  obtained  in  Section  II,  and  these  results  are  compared 
to  the  analytical  results  obtained  in  Section  III  in  the  weak  laser  field 
limit,  i.e.,  e  «  1. 

Equations  (32a, b)  were  numerically  solved  in  order  to  verify  our 

approximation  | AE |  »  e  used  in  deriving  the  constant  of  motion.  Figure  3  is 

a  plot  of  |AEj  as  a  function  of  time  for  the  parameters,  e  m  0.01,  G  *  0  and 

three  different  normalized  frequencies:  f  -  0.96,  f  *  0.98  and  f  *  1.0.  The 

curves  show  that  the  approximation,  | AE }  »  e,  is  indeed  well-satisfied.  The 

curves  for  |AE|  in  Fig.  3  are  periodic  in  time  and  show  that  the  plasma  wave 

periodically  exchanges  its  energy  with  the  laser  field.  The  plasma  build-up 

time  is  longest  when  the  frequency  is  fQ  ,  and  is  29  plasma  periods  with 

f  *  1  and  e  *  0.01,  agreeing  with  the  calculations  in  Section  III.  C. 

Figure  4  is  a  plot  of  the  amplitude  of  E  as  a  function  of  time  for  the 

complete  equations  (18-20),  with  the  same  parameters  as  used  in  Fig.  3.  The 

laser  power  was  built-up  gradually  over  three  periods  of  the  laser  beat 

frequency.  We  used  8  ■  0.9999  for  the  purpose  of  comparing  with  analytical 

& 

results.  In  this  case,  the  amplitudes  of  the  plasma  electric  fields  are  only 
slightly  changed,  showing  that  the  analytical  equations  are  excellent  In  the 
small  e  limit. 

In  Fig.  5,  the  phase  velocities  given  by  the  analytical  results  (solid 

curve)  are  compared  to  the  numerical  values  (dashed  curve)  for  e  *  0.01  and 

f  ■  0.98.  The  shift  of  the  two  curves  is  due  to  the  three  periods  of  laser 

build-up  time  in  the  numerical  calculation.  The  phase  velocity  v  .  is  at  a 

ph 

local  minimum  when  | AE |  is  maximum  and  decreases  as  AE  increases. 

The  plots  of  the  numerically  calculated  peak  amplitude  of  E  (dashed 

1/3 

curve)  and  the  analytical  expression  AE  -  4(e/3)  (solid  curve)  are 

max 
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plotted  as  a  function  of  e  in  Fig.  6.  The  plots  of  the  normalized  laser 

beat  frequency  f  (dashed  curve),  at  which  the  largest  accelerating  electric 

field  is  numerically  obtained  and  the  analytical  expression  for  fQ  (solid 

2/3 

curve)  are  shown  in  Fig.  7  as  a  function  of  e 

The  normalized  peak  amplitudes  of  the  accelerating  electric  field 
for  e  *  0.01,  0.04  and  0.16  in  the  frequency  range  0.8  <  f  <  1.1  are  given  in 
Fig.  8.  The  phase  velocities  associated  with  the  peak  amplitude  are  plotted 
in  Fig.  9.  The  dashed  curves  are  the  numerical  results  obtained  from  the 
complete  nonlinear  and  fully  relativistic  equations  (18-20),  while  solid 
curves  are  the  analytical  results  obtained  from  Eq .  (34).  The  agreement 
for  small  e  is  excellent. 

We  note  that  the  analytical  solutions  for  the  amplitude  of  the  plasma 
wave  and  the  phase  velocity  have  discontinuities  at  fopt*  For  f  >  fopt , 

AE  has  one  real  root  and  the  phase  velocity  associated  with  the  peak  electric 
field  is  less  than  the  speed  of  light.  For  f  <  fQpt ,  AE  has  three  real 
roots.  The  smallest  value  of  the  root  is  closest  to  the  numerical  result. 

For  f  <  fQ  ,  the  phase  velocity  associated  with  the  peak  electric  field  is 
generally  greater  than  the  speed  of  light. 

In  the  region  marked  by  (*■*»  =  »)  in  Fig.  8,  the  analysis  is  not 
applicable  because  of  particle  mixing.  Here,  the  large  amplitude  oscillations 
cause  the  electrons  to  become  relativistic  resulting  in  wave  steepening,  1  ’ ^ 
which  phenomenon  is  illustrated  in  Fig.  10.  Two  curves  of  the  normalized 
electric  field  at  two  Instants  of  time  are  plotted  for  one  wavelength  of  the 
laser  beat  wave  with  parameters  e  -  0.16  and  f  »  0.925,  the  point  is  marked  by 
(*)  on  Fig.  8.  The  curve  with  the  circular  dots  showing  wave  steepening,  is 
the  normalized  electric  field  just  before  wave  breaking.  On  the  other  hand, 
the  wave  profile  is  almost  sinusoidal  when  the  amplitude  of  the  electric  field 


*.  *.  -.  •  *.  %  , 


V  \  \  *. 
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T 


T 


is  small,  i.e.,  | AE J << 1 ,  as  illustrated  by  the  curve  with  the  crosses  in  Fig. 

10.  As  | AE|  becomes  Larger  than  1.1,  the  electric  field  in  our  representation 

becomes  multivalued  and  the  Lagrangian  model  breaks  down.  Since  the  electric 

field  produced  by  a  turbulent  plasma  is  unlikely  to  be  desirable  for 

acceleration  of  electrons,  the  upper  limit  for  |AE|  is  AE  »  1.1.  The 

Lim 

minimum  e  necessary  to  obtain  AE^m  *-s  approximately  e  ^  “  0.06. 

To  illustrate  the  relativistic  phenomenon.  Fig.  11  shows  a  plot  of 
2  ~l/2 

y  »  (l  -  (v  /c)  1  for  one  wavelength  of  the  laser  beat  wave  for  the  same 
z  z 

parameters  and  at  the  same  instants  in  time  as  Fig.  10. 

Numerically  calculated  time  evolution  of  the  electric  field  and  the 
corresponding  phase  velocity  for  f  -  0.925  are  plotted  in  Figs.  12  and  13 
respectively,  where  e  is  adiabatically  increased  from  0  to  0.16  in  three 
plasma  periods. 

Since  the  saturated  oscillatory  electric  field  amplitudes  and  phase 
velocities  are  not  desirable  for  accelerating  electrons,  the  laser  pulse 
duration  should  be  chosen  to  equal  the  plasma  build-up  time.  After  the  laser 
pulses  pass  through  the  plasma,  the  plasma  wave  will  continue  to  oscillate 
until  disrupted  by  various  instabilities.  Figure  14  plots  the  temporal 
profiles  of  the  normalized  stimulated  electric  field  amplitude  (solid  curve) 
for  f  ■  0.925  and  a  short  laser  beat  wave  pulse  e(t)  (dashed  curve).  The 
corresponding  phase  velocity  is  shown  in  Fig.  15.  After  the  laser  pulse 
passed,  the  amplitude  and  the  phase  velocity  of  the  plasma  oscillation 
remained  roughly  constant. 

Next,  we  examined  the  effect  of  the  perpendicular  magnetic  field  on  the 
plasma  oscillation  in  the  surfatron  configuration.  An  imposed  transverse 
magnetic  field  can  increase  the  total  electron  energy  by  maintaining 

22 

synchronism  while  accelerating  the  electrons  in  the  transverse  direction. 
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The  analytical  calculation  shows  that  the  transverse  magnetic  field  has  a 


higher  order  effect  on  the  plasma  dynamics.  The  numerical  result  of  the  peak 
electric  field  and  the  corresponding  phase  velocity  are  plotted  as  a  function 
of  G  for  £  =  0.16  and  f  =  0.925  in  Fig.  16,  and  results  changed  little  for 
0  <  G  <  1.  The  imposed  transverse  magnetic  field  increases  the  effective 
plasma  frequency,  but  has  little  effect  on  the  dynamics  of  the  plasma  wave. 
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V.  Summary 

We  have  obtained  nonlinear,  fully  relativistic  results  for  the  plasma 
waves  excited  by  the  beating  of  two  laser  beams.  We  found  that  the  effective 
plasma  frequency  is  a  function  of  the  laser  power  as  'well  as  the  imposed 
transverse  magnetic  field  in  the  surfatron  configuration.  In  the  ideal 
situation  analyzed  here,  the  amplitude  of  plasma  waves  becomes  oscillatory. 
Since  it  is  desirable  to  maintain  the  accelerating  electric  field  at  the 
largest  value,  the  laser  pulse  duration  should  be  approximately  equal  to  the 
plasma  wave  build-up  time. 

In  the  weak  laser  power  limit,  we  obtained  analytical  results  for  the 

saturated  plasma  wave  for  a  range  of  frequencies  around  the  effective  plasma 

frequency.  As  the  laser  power  increases,  the  maximum  amplitude  of  the  plasma 

1/3 

electric  field  increases  as  e  ,  confirming  the  previous  work  of  Rosenbluth 
and  Liu^  at  exact  resonance  when  the  transverse  motion  is  nonrelativistlc. 


i.e.,  y^0=  1,  and  the  corresponding  phase  velocity  decreases,  i.e.,  v^  - 
2/3 

vph  sca^es  as  E  •  The  maximum  accelerating  field  is  achieved  when  the  laser 
beat  frequency  is  less  than  the  effective  plasma  frequency.  The  difference 

between  the  effective  plasma  frequency  and  the  optimal  laser  beat  frequency  is 

2/3  -2/3 

proportional  to  e  .  The  plasma  build-up  time  is  proportional  to  e 

Given  a  plasma  density,  the  desirable  range  of  normalized  laser  beat 

frequency  for  growth  of  large  amplitude  plasma  waves  is  small,  i.e., 

1.05  >  f  >  f  .  If  the  laser  beat  frequency  is  given,  this  condition  can  be 
-  ~  opt 

translated  to  plasma  density  requirements.  Defining  nQ  -  An  +  n^j,  where 
nf„i  is  the  ambient  plasma  density  that  will  provide  exact  resonance  at  f*l, 

,  /  >  /  2  ,  ^  n2 ,  2  ,1/2  .  ,  /,  i  i 2  ,  .1/2 

i.e.,  a)  (n,  .)  -  (u)  y.  +  0  /y  )  =  Acu,  where  w _ *  (4ir  e  n,  ,/m  ) 

s  f*l  po  io  io  po  f*I  o 

The  density  criteria  for  stimulating  large  amplitude  plasma  waves  is 
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(47) 


-  0.1  (1  +  G  )  <  — —  <  1.1  e2/3  (1  +  G  ), 
o  ~  nf,1  - 

2  2  2  —1 

where  Gq  *  (n  /y^Q) ( t«>  /y^)  *  We  now  consider  two  examples  of  plasma 

density  variation  requirements  without  imposed  transverse  magnetic  field, 
i.e.,  Gq  ■  0.  The  largest  tolerable  plasma  density  variation  is  about  25% 
for  £  ■  0.06.  The  tolerable  density  variation  is  reduced  to  15%  for 
e  *  0.01.  The  application  of  transverse  magnetic  field  Bq  in  the  surfatron 
configuration  can  increase  the  allowable  density  fluctuation  limit  by  a  factor 
(1  +  G  ).  In  an  experimental  situation,  not  only  overcomes  the  problem  of 

desynchronization,  but  allows  more  flexibility  in  the  tuning  of  plasma 
density. 

A  comparison  of  the  numerical  results  from  the  complete  equations  with 
the  simple  analytical  results  for  the  weak  laser  power  limit  is  excellent 
for  e  «  0.1,  and  is  in  fair  agreement  for  larger  e.  Numerical  results  show 
that  the  largest  amplitude  for  the  accelerating  electric  field  in  cgs  units 
without  electrons  overtaking  each  other  is 


to  me 


The  laser  power  required  to  reach  this  value  without  wave  breaking 
is  e  -  0.06. 

For  the  purpose  of  accelerating  electrons,  it  is  desirable  not  only  to 

have  the  largest  accelerating  electric  field  but  also  have  v  ,/c  less  than  but 

on 

close  to  unity.  Instead,  we  find  that  as  the  amplitude  of  plasma  waves  get 
larger,  the  corresponding  phase  velocities  become  smaller.  .An  applied 
transverse  magnetic  field  can  overcome  the  problem  of  desynchronization  of  the 
accelerating  electrons  in  the  accelerative  electric  field.  Since  the 
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transverse  magnetic  field  only  modifies  the  effective  plasma  frequency,  but 
has  little  effect  on  the  plasma  dynamics,  the  surfatron  configuration  may  be 
the  desirable  way  to  operate  the  laser  beat  wave  accelerator. 

Finally,  we  would  like  to  point  out  that  the  laser  plasma  interaction 
contains  a  rich  source  of  instabilities,  many  of  which  may  be  detrimental  for 
the  formation  of  the  large  amplitude  plasma  waves  studied  in  this  paper  for 
the  laser  beat  wave  accelerating  scheme.  Some  of  the  processes^  that  have 
large  growth  rates  are  three-wave  forward  Raman  scattering,  four-wave  forward 
Raman  scattering,  and  processes  associated  with  background  ions.  Other  areas 
requiring  investigations  are  the  effects  of  the  transverse  Weibel  instability 
induced  by  energy  anisotropy,  the  influence  of  kinetic  effects,  self-focusing 
of  laser  radiation  and  f ilamentation.  Detailed  studies  of  them  are  necessary 
to  evaluate  feasibility  of  the  long  term  goals  of  the  laser  beat  wave 
accelerating  scheme. 
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Appendix  A:  Alternative  View  of  Equations  (32a, b) 


The  simplified  equations  in  (32a,b)  can  be  written  in  the  form 


dx  ,  c2  3  ,  2  ^  2,, 

"’77  y  }l  y> 


(  A-l  ) 


dy  e  1  , ,  ,2  3.2  2  . 

m  T - 1  "  f  -  q  (x  +  y  )  1  x. 


dT  4  2  f2 


(A-2) 


where  x  *  AE  cos  9,  y  *  AE  sin  9  and  the  initial  conditions  are  x  ■  y  ■  0  at 
T  ■  0.  The  amplitude  and  the  phase  of  the  electron  displacement  are 


AE  -  (x  +  y  ) 


2  -  ••2'1/2  and  9  -  tan"1  (-2-). 


(A-3a,b) 


Figure  17  is  a  plot  of  x  and  y  for  e  ■  0.01  and  three  different  normalized  laser 

beat  frequencies:  f  -  0.96,  0.98  and  1.0.  For  f  >  fQpt ,  the  enclosed  area  is  to 

the  left  rf  x  -  0,  and  to  the  right  of  x  -  0  for  f  <  f  . 
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Fig.  3  Plot  of  | AE|  as  a  function  of  number  of  laser  beat  wave  periods 

for  e  ■  0.01,  G  ■  0  and  three  different  beat  frequency  parameters 


f  *  0.96,  0.98  and  1.0. 


Amplitude  Of  Electric  Field  (E) 


Fig.  4  Normalized  amplitude  of  the  accelerating  electric  field 

E  as  a  function  of  the  number  of  laser  beat  wave  periods  obtained 
with  the  fully  relativistic  equations  for  laser  powers  built-up  over 
three  laser  periods. 
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Laser  Power  Parameter  (€) 
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Fig.  5  Phase  velocities  obtained  from  the  simplified  equations  (solid  curve) 
and  numerically  obtained  from  the  fully  relativistic  nonlinear 
equations  (dashed  curve)  for  e  -  0.01,  G  ■  0  and  f  -  0.98. 
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Of  Electric  Field 


ig.  6  Comparison  of  based  on  the  analytical  result  (solid  curve),  and 


the  peak  amplitude  of  E  numerically  obtained  from  the  fully 
relativistic  nonlinear  equations  (dashed  curves). 


Fig.  7  Comparison  of  the  frequency  at  which  the  largest  accelerating 
electric  field  is  obtained.  The  solid  curve  is  a  plot  jf  f 

opt 

obtained  analytically  and  the  dashed  curve  is  found  numerically  from 
the  fully  relativistic  nonlinear  equations. 


Peak  Am 


Fig.  8  Normalized  peak  amplitude  of  the  accelerating  electric  field 
for  £  ■  0.01,  0.04  and  0.16  within  the  normalized  laser  beat 
frequency  range  0.8  <  f  <  1.1.  The  solid  curves  are  analytical 
results  of  peak  AE  obtained  from  Eq .  (34)  and  the  dashed  curves  are 
numerical  results  of  peak  amplitude  of  E  obtained  from  the  complete 
analysis  of  Eqs.  (18)-(20). 
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Fig.  9  Plots  of  the  phase  velocity  associated  with  the  peak  electric  field 
for  £  *  0.01,  0.04  and  0.16  within  the  normalized  beat  frequency 
range  0.8  <  f  <  1.1.  The  solid  curves  are  obtained  from  the 


analytical  equations  and  the  dashed  curves  are  obtained  from  the 
complete  equations  at  the  maximum  of  E. 
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Fig.  10  Normalized  electric  field  for  one  period  of  the  laser  beat  wave  at 
two  different  amplitudes  obtained  at  two  different  instants  In  time 
with  e  M  0.16,  f  =*  0.925  and  G  *  0,  the  point  is  marked  by  (*)  in 
Fig.  8. 
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Number  Of  Laser  Beat  Wave  Periods 

Fig.  15  The  phase  velocity  associated  with  the  pulsed  laser  beam  shown  in 
Fig.  14. 
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Transverse  Magnetic  Field  Parameter  (G) 

Fig.  16  Plot  of  the  maximum  amplitude  of  the  normalized  electric  field  E 


(solid  curve)  and  the  corresponding  phase  velocity  variation 
v  , /v  (dashed  curve)  as  a  function  of  G  for  e  *0.16  and  f  *  0.925 
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Fig.  17  Plot  of  AE  sin  9  versus  AE  cos  9  for  e  *  0.01,  and  three  different 
normalized  laser  beat  frequencies. 
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